Coherent control of intersubband optical bistability in quantum wells 
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We present a study of the nonlinear intersubband (ISB) response of conduction electrons in a 
GaAs/Alo.sGao.rAs quantum well to strong THz radiation, using a density-matrix approach com- 
bined with time-dependent density-functional theory. We demonstrate coherent control of ISB 
optical bistability, using THz control pulses to induce picosecond switching between the bistable 
states. The switching speed is determined by the ISB relaxation and decoherence times, Ti and T2. 



Intersubband (ISB) transitions in semiconductor quan- 
tum wells take place on a meV energy scale and are 
therefore attractive for THz device applications.^ Many 
ISB effects of practical interest occur in the nonlinear 
regime, such as second- and third-harmonic generation?^ 
intensity-dependent saturation of photoabsorption?^ di- 
rectional control over photocurrentsj^ generation of ul- 
trashort THz pulses,^ plasma instability,^ or optical 
bistability.*'^ Inspired by the photoabsorption experi- 
ments by Craig et al.^ this Letter presents a theoretical 
study of the optical bistability region in a strongly driven, 
modulation n-doped GaAs/Alo.sGao.rAs quantum well. 
We here demonstrate, for the first time to our knowledge, 
that ISB bistability can be manipulated on a picosecond 
time scale by short THz control pulses. This opens up 
new opportunities for experimental study of optical bista- 
bility, which in turn may lead to new THz applications 
such as high-speed all-optical modulators and switches. 

Most previous theoretical studies of nonlinear ISB dy- 
namics were based on density-matrix approaches, ei- 
ther in Hartree approximationiSiiiii^ or using exchange- 
only semiconductor Bloch equationsr^^"* In this work, 
we use a density-matrix approach combined with time- 
dependent density-functional theory, which has the ad- 
vantage of formal and computational simplicity, while 
including exchange-correlation (xc) many-body effects. 

We describe the conduction subbands in effective-mass 
approximation, where m* = 0.067to and e* = e/-\/e, 
e — 13, are the effective mass and charge for GaAs. 
Initially, the conduction electrons are in the ground 
state, characterized by single-particle states of the form 



^(0 



r) — A ^"e*'iii""iii/;P(2;), with r|| and q|| in the x — y 

plane. The envelope function for the jth subband ipj{z) 
follows self-consistently from a one-dimensional Kohn- 
Sham (KS) equation^i^ with the ground-state density 



n{z) 



(1) 



Here, Ejq^^ = Sj +h q^^/2m*, and Sj and ep are the sub- 
band and Fermi energy levels. We choose the electronic 
sheet density Ng such that only the lowest subband is 
occupied, in which case ep — nTi^ Ng/m* + ei. 

Under the influence of THz driving fields, linearly po- 
larized along z, the time-dependent states have the form 
^fjqii(r,t) = e*iii'"ii'0j(z,i), with initial condition 



In the absence of disorder and 



,qi|(r,io) = *°q||(r). 
phonons, the time evolution of the envelope functions 
follows from a time-dependent KS equation: 
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2m* dz^ 



+ Vdr{z,t) + 

'^conf 

{z)+Vii{z,t)+V^ 



(3) 

Here, Vdriz,t) = e£ z f (t) sm{ujt) describes the driving 
field, with electric field amplitude £, frequency lu, and 
envelope f(t). Wconf (-z) is the confining square-well poten- 
tial, the Hartree potential follows from cfvB_{z,t)/dz^ — 
— 47re*^n(z, t), and we use the time-dependent local- 
density approximation (TDLDA) for the xc potential: 
v^c{z,t) = [dexc("-)/rf?i]fi=n(^,t) (Sxc is the homogeneous 
electron gas xc energy density). The time-dependent den- 
sity n{z,t) follows by substituting tpj{z,t) in Eq. (Q). 

To account for disorder or phonon scattering, Eq. Q 
is replaced by a density-matrix approach. We expand the 
first conduction subband as tpi{z,t) = J2k=i'^k{t)ip'l{z). 
The associated Ni, x A^;, density matrix p has elements 
Pkiit) = cl{t)ci(t) and initial condition pki{to) = Suhk- 
The time evolution of p follows from 



in^ = [h{t),p{t)]-R, 



(4) 



with the relaxation matrix Rki = h[pki{t) — Pki{U))\/Tki- 
For simplicity, Tki = TiSu +T2{1- Ski), where Ti and T2 
are phenomenological relaxation and decoherence times. 

We consider a 40 nm GaAs/Alo.aGao.rAs square quan- 
tum well with Ns = 6.4 x 10^" cm~^, following Craig et 
al^ The lowest subband spacing is £2 — £1 = 8.72 meV, 
and the ISB plasmon frequency is found from linear- 
response theorjii^ as Ticuisb — 9.91 meV. The system has 
9 bound levels, and we take the lowest five to construct 
the density matrix {Nb — 5). We use the ISB scattering 
times Ti = 40 ps and T2 = 3.1 ps, consistent with recent 
measurements of Ti and T2 for similar systems^Siii 

To describe ISB photoabsorption, we propagate Eq. 
Q in the presence of THz driving fields, switched on at 
to over a 5-cycle linear ramp and then kept at constant 
intensity for several hundred ps. From the dipole mo- 
ment d{t) = NsJ^ki' Pki{t) J dz il;'l{z)zip'j'{z) we obtain 
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FIG. 1: ISB photoabsorption and occupation probabil- 
ity p2 of the second subband for a 40 nm GaAs/Alo.aGao.yAs 
quantum well with electron density 6.4 x 10^*^ cm~^ , driven 
by THz fields with intensities as indicated. Lines: two-level 
RWA^ Points: density-matrix calculation [Eq. Q]. 

the photoabsorption cross section (the dissipated power) 
~ w cos{u)t')d{t') dt', where T is one cycle of 
the driving field Vdr- Following the switching on of the 
THz field, <7{uj) fluctuates considerably from one cycle to 
the next, but eventually approaches a stable value as the 
transients settle down. 

Fig. n shows cr(w) calculated with our density-matrix 
formalism, in good agreement with a two-level rotating- 
wave approximation (RWA) modeli£ including xc. For 
low intensities, (t{lo) is a Lorentzian with maximum at 
cjisB and width 2h/T2. Approaching the saturation in- 
tensity lo = cei/^fi^/STrTiTse^KV'glzlOP, the photoab- 
sorption peak shifts to lower energies and changes shape 
(here, Iq = 26.8 W/cm^). The physical reason for this ef- 
fect is that the depolarization shift A = hujisB — (£2 — £1) 
is proportional to the population difference pi — P2- Ini- 
tially, pi = 1 and P2 — 0- Strong driving fields lead to a 
decrease of pi — p2 because of transitions into the second 
and, eventually, higher levels, see the lower panel of Fig. 
n Since this population transfer is most efficient around 
the ISB resonance, the peak of cr(w) shifts more than the 
tails, leading to an asymmetric line shape. 

For intensities greater than 16 W/cm^ we discover re- 
gions of optical bistability in cr(uj): the system has two 
different ways of responding to one and the same driving 
field (the middle branch is unstable). The two bistable 
states are characterized by different amplitudes of the 
density oscillations, and different level populations. 

Whether the system will be in the upper or lower 
bistable state depends on its history. A hysteresis-like 
behavior is observed at fixed intensity and under adia- 
batically slow frequency changes of the driving field: 
Entering the bistability region from the low-frequency 
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FIG. 2: Intensity and length of THz control pulses used 
to switch from the lower to the upper ISB bistable state at 
I = 30 W/cm^ and uj = 9.3 meV driving field (0.445 ps per 
cycle). Contour labels indicate the switching times in ps. No 
switching occurs in the hatched region. 



side, the system continues on the lower branch and then 
jumps up at the end of the bistability region. Entering 
from the other side, the system follows the upper branch. 
The required continuous, adiabatic frequency tuning of 
THz driving fields is difficult to realize in experiment, and 
of little practical use in exploring ISB optical bistability 
for potential applications. Instead, it would be desirable 
to switch rapidly between the two bistable states. 

In the following, we demonstrate coherent control of 
ISB optical bistability by short THz control pulses. This 
method is both rapid and robust, and lends itself for 
experimental implementation. Fig. suggests that a 
switch from the lower to the upper bistable state requires 
a transfer of population into the second subband level. 
The necessary energy can be rapidly injected into the sys- 
tem by a short THz pulse. On the other hand, switching 
down from the upper to the lower bistable state requires 
upper-level population to relax. We therefore expect this 
process to be inherently slower than the up-switch. 

We consider a driving field with / = 30 W/cm^ and 
hu! — 9.3 meV (0.445 ps per cycle) in the middle of 
the bistable region, such that the system is in the lower 
bistable state. At ti, we apply a short pulse with the 
same frequency and in phase with the driving field, and 
with a trapezoidal pulse envelope: linearly turned on and 
off over 5 cycles, constant in between over Nc cycles (the 
precise pulse shape is not essential). We calculate d{t) 
and 17(0;) as above, using our density-matrix formalism. 
After transients and other disturbances induced by the 
pulse have subsided, the system either slowly returns to 
the lower state, or converges towards the upper bistable 
state, in which case we define a "switching time" as that 
time after ti when cr(u;) is converged to within 5% of its 
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FIG. 3: Same as Fig. H but now for switching from the upper 
to the lower ISB bistable state. The THz control pulses are 
phase shifted by tt with respect to the driving field. 

final value. To study the reverse process, we prepare the 
system in the upper bistable state and then apply simi- 
lar control pulses, phase shifted by tt with respect to the 
driving field. After the transients have settled, the sys- 
tem either remains in the upper state, or goes down to 
the lower state. 

Control pulses of various intensity and length Nc were 
tested to determine the conditions for successful up- or 
down-switching, and to find those pulses that induce the 
fastest switching. The results are summarized in Figs. 
13 and |3| The hatched areas indicate the unsuccessful 
pulses: in the case of up-switching, they do not have suf- 
ficient energy to promote enough electrons to the upper 



level, and in the case of down-switching, they are too 
short to allow enough electrons to relax down. 

The limiting curve in Fig. |21 separating the successful 
from the unsuccessful pulses corresponds to pulses of en- 
ergy 1.87 X 10~^° J/cm^. The shortest switching times 
of 15-20 ps are achieved for somewhat higher pulse en- 
ergy, around 3.3 x 10"^" J/cm^. No systematic effort was 
made to further optimize the control pulses. A physical 
limit for switching speed is set by the decoherence time 
T2 (here 3.1 ps), which governs the decay of transients. 

Down-switch from the upper to the lower ISB bistable 
state requires longer control pulses, at least 30 cycles. 
The associated switching times are limited by the ISB re- 
laxation time Ti (here 40 ps). We find that the shortest 
down-switching times of 50-60 ps are achieved by con- 
trol pulses with intensity not far from the driving field, 
30 W/cm^. Due to the phase difference tt, these control 
pulses effectively suppress the driving field, giving elec- 
trons time to relax from the upper level. 

In summary, we have simulated coherent control of ISB 
optical bistability in quantum wells, using short THz con- 
trol pulses to switch rapidly between the bistable states. 
Switching times are in principle limited only by the in- 
trinsic relaxation and decoherence times. Interestingly, a 
shorter relaxation time Ti (i.e., more scattering) implies 
faster switching. More realistic simulations might include 
microscopic theories for Ti and T2 and their intensity de- 
pendence, finite temperatures, and a detailed modeling 
of the absorption profile and waveguide geometry of the 
device. Ultimately, ISB optical bistability may pave the 
way towards new THz active optical elements. 
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